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The high-mass microquasar Cygnus X-3 has been recently detected in a flaring state by the 
gamma-ray satellites Fermi and Agile. In the present contribution, we study the high-energy 
emission from Cygnus X-3 through a model based on the interaction of clumps from the Wolf- 
Rayet wind with the jet. The clumps inside the jet act as obstacles in which shocks are formed 
leading to particle acceleration and non-thermal emission. We model the high energy emission 
produced by the interaction of one clump with the jet and briefly discus the possibility of many 
clumps interacting with the jet. From the characteristics of the considered scenario, the produced 
emission could be flare-like due to discontinuous clump penetration, with the GeV long-term 
activity explained by changes in the wind properties. 
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Figure 1: Left: Sketch of a HMMQ. Right: Sketch of the interaction of one clump with the jet. 



1. Introduction 

The mass loss in massive stars is thought to take place via supersonic inhomogeneous winds. 
Considerable observational evidence supports the idea that the wind structure is clumpy [P, al- 
though the properties of clumps are not well-known as a consequence of the very high spatial 
resolution necessary for clump detection. Some massive stars are accompanied by a compact ob- 
ject and present transfer of matter to the latter, forming an accretion disc. When bipolar relativistic 
outflows are formed, these sources are called high-mass microquasars (HMMQs). In Figure [I] (left), 
a sketch of a HMMQ is shown. 

Cygnus X-3 is a HMMQ composed by a Wolf-Rayet (WR) star and a compact object, both 
being separated a distance a ~ 3 x 10 11 cm. The nature of the compact object has not been deter- 
mined yet, but recent studies suggest that it may be a black hole (e.g. Located at a distance 
~ 7 — 9 kpc, the system has been recently detected in a recurrent flaring state at high-energy (HE) 
gamma rays [^J ffl. The observed emission spectrum is a power- law with index ~ 1 .7 and luminos- 
ity > 10 36 erg s _1 . The typical duration of the flares was ~ 2 days PJ]. 

In this contribution, we present a model to explain the HE emission detected in Cygnus X-3 
based on the interaction between its jets and wind inhomogeneities [Q]. Some clumps can even- 
tually penetrate in the jet, leading to transient non-thermal activity that may release a significant 
fraction of the jet kinetic luminosity, Lj, in the form of synchrotron, inverse Compton (IC), and 
proton-proton (pp) emission. 

2. Jet-clump interaction in Cygnus X-3 

In order to study the physical processes of the jet-clump interaction, we consider the collision 
of a single clump with the jet, as is shown in Figure [I] (right). A huge density contrast % = p c / pj 
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(where p c and pj are the clump and jet densities, respectively) allows the clump to fully cross the 
boundary of the jet and penetrate into it. In the context of this work, thermal conduction, magnetic 
fields and gravitational forces are not dynamically relevant for the jet-clump interaction and have 
been neglected. 

We consider that a clump with size R® = 10 10 cm (< 0.1/?*) reaches one of the HMMQ jets 
moving at the terminal wind velocity = 2500 km s~ 1 . At a distance r = Vz 2 + a 2 from the star, 
the clump density can be estimated through p c = p w /f, where / < 1 is the filling factor of the 
clumpy wind, and p w ~ M / (4nr 2 Vco) is the wind density, where M~ 10 _5 M Q yr _1 is the typical 
mass loss rate of WR stars. We assume that the jet has a velocity vj ~ c/3 (i.e. Lorentz factor 
r ~ 1.06) and a kinetic luminosity L s = 10 38 erg s . We fix the relation between the height (z) 
and the width (Rj) of the jet as Rj = O.lz. 

The penetration time of the clump into the jet is determined by t c ~ 2R° c /v c ~ 100 s, where the 
clump velocity is v c = Voo. As a consequence of the interaction of the jet material with the clump, 
two shocks form. One of these shocks propagates back in the jet with a velocity Vb s ~ vj, forming 
a bow shock (see Figure [j], right). This bow shock reaches the steady state configuration in a time 
?bs ~ ZbsAbs, w * tn tne sta g na ti° n point located at a distance Zb s ~ R®/5 from the clump. On the 
other hand, a shock propagates in the clump at a velocity v sc ~ vj/ 1 y/% and, in a time ? sc ~ 2R c /v sc , 
the whole clump is shocked. During the shock passage, the clump is compressed in the direction 
of the shock velocity, z, and when the whole clump is shocked, its size in the z-direction is ~ 
However, in t > ? sc the clump starts to expand at the sound velocity of the shocked material. At 
the same time, the acceleration exerted by the jet to the clump increases with R c . Assuming that 
the expansion is spherical, and neglecting the acceleration exerted by the jet on the clump, we 
obtain the following expression for the clump size: R c (t) ~ /?^/(l — 0.4f/? sc ) 2 [n], which gives a 
characteristic clump expansion time 

^exp ~ 3.5? sc . In Figure ^ (left), the evolution of R c with time 

is shown. 

As a consequence of the acceleration exerted by the jet on the clump, Rayleigh-Taylor (RT) 
instabilities can develop in the contact discontinuity. Kelvin-Helmholtz (KH) instability can also 
grow as a result of the high relative velocity between the jet shocked material and the clump. 
Numerical simulations show that RT and KH instabilities grow sufficiently to destroy the clump 
(i.e. up to a wavelength ~ R c ) in a timescale £rt ~ ?kh ~ 5? sc [|7|]. If instabilities do not grow 
fast enough, and clump expansion and acceleration are too slow, the permanence of the clump 
into the jet will be determined by the passage time of the clump through the jet: £j ~ 2R-Jv c . In 
Figure || (right), the main dynamical timescales are plotted. 

2.1 Interaction height 

Considering the previous analysis on the dynamics of the jet-clump interaction, one can deter- 
mine the minimum interaction height, z™J n , at which one clump can completely penetrate into the 
jet. Being the lifetime of clumps into the jet fijf e > ? sc , to determine the interaction height z; nt we 
impose that ? sc > t c . This constraint is fulfilled if z > AVz 2 + a 2 , where 

— 1/2 

0.5 f — ) f -L ) f ._,„ J , ) f (io-5M^/yr) (2500km/s) 

(2.1) 



0.01 J \c/3 J V 10 38 erg/s 
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If, for a certain election of clump and jet parameters, A results larger than 1, the constrain f sc > t c 
is no satisfied for any value of z and clumps will be destroyed by the jet before fully penetrating 
into it. However, for the parameters chosen in this work, A ~ 0.5 and in such a case, the minimum 
interaction height results z|™ n ~ Aa/\J\ —A 2 ~ a/2. 

3. Non-thermal processes 

To estimate the non-thermal emission produced by the interaction of one clump with the jet, we 
fix the interaction height at zm = a ~ 2z™J". At zi nt , the density of the clump is p c ~ 10~ 12 gr cm~ 3 
(i.e. the number density is n c = p c / (4m p ) ~2x 10 11 cm~ 3 , where m p is the proton mass). 

3.1 Particle acceleration 

In the bow shock, we assume that particles are accelerated up to relativistic energies being 
injected in the downstream region following a distribution Q e p oc E~ p (e and p for electrons and 
protons, respectively) and with a luminosity L nt = T] nt (/? c //?j) 2 Lj, where we fix T7 nt = 0.1. For the 
parameters assumed in this work, 

*- ~ 10 " (£) (lO^) 2 G = 3 x^O" en, )' 2 (lO^) ^ K <3 » 

In addition, considering that the magnetic energy density in the bow-shock region, Ub, is a fraction 
tjb of the non-thermal one, U at , we obtain a magnetic field B ~ y^rjB 2 x 10 3 G. (Both Tj nt and tjb 
are free parameters in our work.) 

The main radiative losses that affect the evolution of the non-thermal particles are synchrotron 
radiation, and synchrotron self-Compton (SSC) and external Compton (EC) scattering. For EC 
cooling we have considered target photons provided by the WR star, that has a luminosity ~ 
10 39 erg s _1 and a temperature ~ 10 5 K. The energy density of stellar photons at zj nt is t/* ~ 
3 x 10 4 erg cm~ 3 . At electron energies E e > 25 GeV, EC scattering with stellar photons of energy 
£ph* ~ 3AT* (where K is the Boltzmann constant) occurs in the Klein-Nishina regime. In addition 
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Figure 3: Left: Timescales of the acceleration and lepton cooling losses (synchrotron, SSC, EC and relativis- 
tic Bremsstrahlung). Diffusion and advection escape losses are also shown. Right: Spectral energy distribu- 
tion produced by one clump of = 10 10 cm interacting with the jet at Zmt ~ a- At energies £ p h > 10 GeV, 
the emission is strongly absorbed. Both (left and right) plots have been produced taking B~2x 10 2 G. 

to radiative losses, electrons are advected away from the emitter on a time ~ 4R c /v y In Figure || 
(left), the electron cooling timescales are plotted together with the acceleration timescale. As seen 
in the figure, for T7b = 10~ 2 the maximum energy is determined by synchrotron losses, yielding 
E™ ax ~ 1 TeV. To obtain the distribution N e of relativistic electrons we solve the kinetic equation [||] 
for a one-zone emitter, taking into account the losses mentioned above, and also the synchrotron 
photons (as targets for SSC). The steady state is reached on a time <C t sc . N e (E e ) has a break at the 
energy ~ 0.1 GeV due to electron advection escape. 

Protons can lose energy via pp interactions in the bow-shock region, but diffusion losses dom- 
inate, constraining the maximum energy to E™ ax ~ 2.5 x lO 2 ^^ TeV (for T7 nt = 0.1). The most 
energetic protons, E p > 0. l£"™ ax , can diffuse up to the clump before escaping through advection, 
producing gamma rays via pp interactions. 

3.2 High-energy emission 

In the bow-shock region, using standard formulas [Q| and N e , we estimate the synchrotron, 
SSC and EC emission; relativistic Bremmstrahlung and pp interactions are negligible there. In 
the clump, if energetic protons that arrive from the bow shock are not confined, they will escape 
from the clump in a time t c \ ~ Rq/c before radiating a significant part of their energy. Considering 
that the energy distribution of relativistic protons in the clump is N p ~ Q p t c b we estimate the pp 
emission following the formulae given in JToj] . In Figure [| (right), we show the spectral energy 
distribution (SED) composed by the most important radiative processes in the bow-shock region 
(synchrotron, SSC and EC) and in the clump (pp). The synchrotron emission is self-absorbed at 
energies E p h < 10~ 4 eV, and at photon energies E p h > 10 GeV absorption due to electron-positron 
pair creation in the star photon field is relevant. The achieved luminosity at energies ~ 0. 1 — 10 GeV 
is Lec ~2x 10 35 erg s _1 . Regarding to the pp emission, the estimate obtained in this work is a 
lower-limit. If the protons were confined in the clump, the radiated emission would be larger. 

The SED shown in Figure ^ (right) has been calculated neglecting the expansion of the clump 
inside the jet (see Figure ^, right). If we consider that after ? sc the cloud expands up to a size much 
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larger than = 10 cm, the available non-thermal luminosity will be significantly larger and, as a 
consequence, the synchrotron and IC luminosity will be higher. On the other hand, the pp emission 
produced in the expanded clump will be lower, because the clump density decreases faster than the 
increase of the clump surface: n c /n® °< (R®/R c ) 3 . 

Finally, note that in Figure |3] (right) we show the result of the interaction of only one clump 
with the jet, but many clumps could simultaneously interact with the jet. The non-thermal radiation 
fluxes grow linearly with the number of clumps interacting with the jet, N c j, as long as too many 
clumps do not disrupt the jet. The corresponding SED will be also higher. 

4. Discussion 

The total luminosity produced by jet-clump interactions depends on the number of clumps 
inside the jet, and on the size and density of each one. The number of clumps in the jet can 
be estimated as N c j ~ /jVj/V c , where ft < f is the filling factor of clumps into the jet, and Vj 
and V c are the jet and the clump volume, respectively. Considering Vj up to z ~ 2 a, V c ~ nR% 3 
(where R® = 10 10 cm) and f\= f = 0.01, we obtain N c j ~ 0.5 in Cygnus X-3. However, the real 
number of clumps into the jet will be smaller than the previous estimate, because f s < f as a 
consequence of the clump destruction during the penetration process and inside the jet due to RT 
and KH instabilities. If N$ < 1, the emission produced by the jet-clump interaction will be like a 
flare with a duration determined by the lifetime of the clumps into the jet. 

On the other hand, if we consider clumps with a radius R® = 10 9 cm (< 0.01 /?*), the number 
of clumps into the jet grows up to ~ 10 4 . Then, the simultaneous interaction of ~ 10 4 clumps with 
the jet will produce significantly more luminosity than the emission produced by the interaction of 
only one clump, that is ~ 10 34 erg s _1 at Zmt ~ «• In that case, the final spectrum will be persistent, 
with a flickering due to the continuous interactions of many small clumps. 

However, the jet can be disrupted by the too small many clumps, in particular if the section of 
clumps interacting within Az is comparable with the section of the jet: LR^ ~ /?? . This will be the 
case, neglecting clump expansion, for /?[? < 10 9 cm. Note that periods of smaller clump size may 
lead to jet disruption. Also, changes in the wind could lead to more clumps inside the jet and more 
effective tapping of the jet luminosity. 
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